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ABSTRACT 

Wc consider the energy-driven stochastic state vector reduction equa- 
tion for the density matrix, which for pure state density matrices can be writ- 
ten in two equivalent forms. We use these forms to discuss the decouphng of 
the noise terms for independent subsystems, and to construct "environmen- 
tal" stochastic density matrices whose time-independent expectations are 
the usual quantum statistical distributions. We then consider a measure- 
ment apparatus weakly coupled to an external environment, and show that 
in mean field (Hartree) approximation the stochastic equation separates into 
independent equations for the apparatus and environment, with the Hamil- 
tonian for the apparatus augmented by the environmental expectation of 
the interaction Hamiltonian. We use the Hartree approximated equation to 
study a simple accretion model for the interaction of the apparatus with its 
environment, as part of a more general discussion of when the stochastic dy- 
namics predicts state vector reduction, and when it predicts the maintenance 
of coherence. We also discuss the magnitude of decoherence effects acting 
during the reduction process. Our analysis supports the suggestion that a 
measurement takes place when the different outcomes are characterized by 
sufficiently distinct environmental interactions for the reduction process to 
be rapidly driven to completion. 
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1. Introduction 

Understanding the measurement process has been a persistent problem since the 
inception of quantum mechanics. In the orthodox Copenhagen interpretation, measurements 
are accounted for by invoking a layer of classical, non-quantum mechanical reality; attempts 
to extend quantum mechanics to include the measuring apparatus itself lead to quandaries, 
such as the famous Schrodinger cat paradox. One approach to this problem that has been 
much studied recently [1-8] postulates that the Schrodinger equation is only an approximate 
description of reality, and must be modified by small, nonlinear stochastic terms. These terms 
drive the state vector reduction process, and account for the non-observation of macroscopic 
quantum superpositions in measurement situations. 

The proposal that a stochastic, nonlinear Schrodinger equation provides the phe- 
nomenology of quantum measurement passes a number of consistency tests. In its energy- 
driven form, it leads exactly to probabilities given by the Born rule [4,8,9], and for mea- 
surements on degenerate systems leads to the Liiders projection rule [9]. There are plau- 
sible arguments [7,8], to be elaborated on here, that with a Planckian magnitude of the 
stochastic term, coherence is maintained where observed experimentally, while state vector 
reduction proceeds for measurement situations where discrete outcomes are observed. How- 
ever, the stochastic Schrodinger equation is inherently nonrelativistic [10,11], involving the 
same stochastic differential at all spatial locations. This raises the issue of whether it is 
consistent with clustering - put simply, does the reduction of the state vector in a localized 
measuring apparatus proceed independently of what goes on far away from the laboratory? 
An affirmative answer to this question was given [8] under the assumption that all of the 
universe is governed by the pure state stochastic reduction equation. In this paper we extend 
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this analysis in a number of directions, with the aim of understanding in greater detail the 
stochastic evolution of a "measurement" system coupled to its environment. 

Our discussion is organized as follows. In Sec. 2 we give two equivalent forms of 
the Ito noise term in the stochastic evolution equation for a pure state density matrix, and 
use these to discuss clustering for disjoint subsystems. In Sec. 3 we use one of these forms 
to prove the existence of pure state density matrices whose stochastic expectation gives the 
standard quantum statistical distributions. We also give a mixed state generalization of 
these results that is relevant when the Hamiltonian has degeneracies. In Sec. 4 we consider 
a "measurement" subsystem weakly coupled to an "environment" subsystem, obeying overall 
the density matrix stochastic evolution equation, and derive the corresponding single system 
Hartree or mean field stochastic equations for the measurement and environment subsystems, 
working to first order accuracy in the interaction Hamiltonian. 

In Sec. 5 we give a survey of under what circumstances the stochastic evolution 
equation predicts state vector reduction, based on Planckian estimates for the magnitude of 
the stochastic term. We show that when the energy spread between superimposed states is 
small, coherence is maintained, in agreement with recent experiments on quantum coherence 
in large systems. On the other hand, when the energy fluctuations in a measurement sys- 
tem are large enough, state vector reduction proceeds rapidly to completion. We consider 
three types of energy fluctuations: thermal fluctuations, shot effect fluctuations in electric 
currents, and surface accretion fluctuations. For the latter, we use the mean fleld approx- 
imation of Sec. 4 to construct a simple model for accretion processes, which motivates a 
quantitative discussion of their influence on the measurement process in both terrestrial and 
extraterrestrial environments. In Sec. 6 we discuss the coexistence of standard decoherence 
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mechanisms with the stochastic reduction process. In Sec. 7 we state our conclusions regard- 
ing the imphcations of this analysis for the measurement process in quantum mechanics. In 
an Appendix, we discuss a coherent state variant of the accretion model of Sec. 5, in which 
reduction can proceed to coherent states. 

2. Stochastic Density Matrix Equations and Clustering 

We begin by recalling some formulas from the theory of stochastic Schrodinger equa- 
tions [1-8]. Letting |x) be a normalized state vector, the standard stochastic evolution 

("quantum state diffusion") equation for \x) takes the form 

d\x) = [adt + (5dWt]\x) , (la) 
with dWt a real Ito stochastic differential obeying 

dW^ = dt , dWtdt = , (16) 

and with 

a = -tH-la'[A-{A)r , 

. ' (Ic) 

f3=-a[A-{A)] , 

where cr is a numerical parameter governing the strength of the stochastic and drift terms, 
and A is a. self-adjoint operator with expectation {A) in the state 

{A) = {x\A\x) . (Id) 

The energy-driven case of the stochastic evolution is obtained by taking A to be the Hamil- 
tonian H, which we shall assume henceforth. It is straightforward to show that the evolution 
of Eqs.(la)-(ld) preserves the normalization of the state vector |x). 



Defining the pure state density matrix p — |x)(x|) it is easy to show that the state 
vector evolution of Eqs. (la)-(ld) imphes that the density matrix evolution is given by 



Direct calculation from Eqs. (la)-(ld) gives the coefficient N{p^ H) of the Ito noise term 
dWt in Eq. (2) in the form 



which by use of the pure state condition (? — p can be written in the equivalent form 



Both of these forms have the property that p^ — p implies that {p, dp} + {dpY — dp, which 
can be rewritten as (p + dp)"^ = P + dp, and so they preserve the pure state condition. 

Let us now consider a system for which the Hamiltonian H is the sum of two Hamilto- 
nians Hi , H2 which depend on disjoint sets of variables, and investigate the conditions under 
which Eqs. (2) and (3a), (3b) admit factorized solutions p — piP2, with pi 2 obeying equa- 
tions of similar form driven by the respective Hamiltonians i?i,2- Substituting H = Hi + H2 
and p — P1P2 into Eqs. (3a), (3b), and using the facts that all variables in set 1 commute 
with all variables in set 2, and that Tr = TriTr2, we find respectively from Eqs. (3a) and 
(3b) that 



N{pip2, Hi + H2) = P2KP1, Hi} - 2piTr2P2Tripiifi] + pi[{p2, H2} - 2p2TripiTr2P2-f^2] 




(2) 



N{p,H)^{p,H}-2pTrpH 



(3a) 




(36) 



(4a) 



N{pip2, Hi + H2) = pI\pu [Pi, ^1]] + P?[P2, [P2, H2]] 



(46) 
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Clustering requires that 

A^(piP2, + H2) = P2Ni{pi, Hi) + piN2{p2, H2) , (5) 

with Ni^2 the restrictions of to the 1,2 subspaces. We see that Eq. (4a) obeys the clustering 
property by virtue of the trace conditions Tripi = 1, Tr2P2 = 1, while Eq. (4b) satisfies the 
clustering property by virtue of the pure state conditions pi = Pi, P2 = P2- 

Let us now examine the clustering properties of the remaining terms in Eq. (2). For 
the left hand side, we find by use of the Ito extension of the chain rule that 

d{p\p2) = p2dp\ + pidp2 + dpidp2 . (6a) 

Thus, in order to have dpi and dp2 obeying equations of the same form as dp but restricted 
to the 1, 2 subspaces, the left hand side should take the form, using Eqs. (lb) and (2), 

d{piP2) = P2dpi + Pldp2 + '^(T^Ni{pi, Hi)N2{p2: H2)dt . (66) 

For the dt terms on the right hand side of Eq. (2), we have 

-i[Hi + H2, piP2]dt-la^[Hi + H2, [Hi + H2, piP2]]dt 

o 

=P2{-z[Hi, pi]dt - la^[Hi, [Hi,pi]]dt} 

I (6c) 

+Pl{-t[H2,P2]dt - -CT^[H2, [H2,P2]]dt} 
o 

~a^[Hi, pi][H2, P2]dt . 
Assuming the conditions for the clustering property of Eq. (5) to hold for the Ito noise term, 
comparing Eqs. (6a)- (6c) we see that the complete density matrix evolution equation will 
cluster if and only if 

Ni(pi,Hi)N2(p2,H2)^-[Hi,pi][H2,p2] . (7) 
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This condition does not hold as in identity for either of the two possible forms for N{p^ H) 
given in Eqs. (3a), (3b), and so the a'^dt or drift term in the stochastic evolution equation 
does couple disjoint systems. 

However, there are two important special cases in which disjoint systems decouple 
asymptotically. The first of these cases corresponds [8] to taking N{p, H) as in Eq. (3b), so 
that Eq. (7) becomes 

[pi, [pi, i/i]][p2, [P2, H^W = -[Hu Pi][H2. P2] . (8a) 

This equation is satisfied, by virtue of both the left and right hand sides vanishing, whenever 
either [pi. Hi] = or [p2, -^^2] = 0, conditions that are respectively obeyed when system 1 or 
system 2 is at the endpoint of the state vector reduction process. In particular, if system 1 
represents a measurement process, and system 2 represents a pure state environment at the 
endpoint of its reduction process, then the stochastic dynamics of system 1 is completely 
independent of the dynamics of its environment. 

A more general case in which disjoint systems decouple asymptotically corresponds 
to taking N{p,H) as in Eq. (3a), but not assuming the pure state condition so that this 
cannot be transformed to Eq. (3b). Equation (7) now becomes 

[{pi. Hi] - 2piTripii/i][{p2, H2} - 2P2TY2P2H2] = -[Hi,pi][H2, P2] ■ {Sb) 

This equation is satisfied, again by virtue of both the left and right hand sides vanishing, 
whenever either pi is a linear combination of projectors on a degenerate submanifold of Hi, 
or p2 is a linear combination of projectors on a degenerate submanifold of H2. For example, 
in the latter case we would have p2i^2 = -f^2P2 = -E2P2 for some degenerate submanifold 
energy E2, together with Tr2P2 = 1, which imply the simultaneous vanishing of {p2, H2} — 
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2p2Tr2P2-f^2 and of [i?2,P2]- Thus, if one were to adopt Eqs. (2) and (3a) as a generalization 
of the density matrix evolution equation to the case of non-pure state density matrices, a 
pure state measurement process decouples from a mixed state environment whenever the 
density matrix for this environment is a linear combination of projectors on a degenerate 
submanifold of its Hamiltonian. 

3. Mctrtingale Construction of the Standard 
Quantum Statistical Distributions 

In order for the measurement system to decouple from its environment, we have seen 
that the environment must be described either by a pure state density matrix that commutes 
with the environment Hamiltonian, or by a mixed state density matrix that is a linear 
combination of projectors on a degenerate submanifold of the environment Hamiltonian 
(with the second case equivalent to the first for a one dimensional submanifold). This raises 
the question of how such a description can be compatible with the usual description of 
equilibrium environments in terms of the standard quantum statistical distributions, which 
are mixed state density matrices p obeying the trace condition Trp = 1, but which do not 
obey either the pure state condition = p or the more general condition that the density 
matrix be a linear combination of projectors on a degenerate Hamiltonian submanifold. The 
answer is that in the theory of stochastic state vector reduction, the role of the usual mixed 
state density matrix is played [7] by the stochastic expectation E[p] and not by p itself. 
Thus an equilibrium environment can be described by a stochastic density matrix that is a 
linear combination of projectors on a degenerate Hamiltonian submanifold, the stochastic 
expectation of which has the form E[p\ — f{H), with / one of the standard quantum 
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statistical distribution functions of the Hamiltonian. Since Eq. (2) implies that E[p\ obeys 
the time evolution equation 

dE[p] = -i[H, E[p]]dt - ^a'[H, [H, E[p]]]dt , (9) 

any E[p\ of the form f{H) is time independent, as expected of the quantum statistical 
distributions. 

To show that there are pure state density matrices with the required expectation, 
we proceed constructively by use of the density matrix evolution equation in the form 

dp = -i[H, p]dt - la^[H, [H, p]]dt + }-a[{p, H} - 2pT^pH]dWt . (10) 
o Z 

Although we derived this equation in Sec. 2 for pure state density matrices, we shall now use 
it, as suggested in the discussion associated with Eq. (8b), as a stochastic evolution equation 
for density matrices p that do not obey the pure state condition. Taking the initial p at 
time t = as po = f{H), we see from Eq. (9), which follows by taking the expectation of 
Eq. (10), that E[p] — f{H) for all times. Also, since Eq. (10) only involves the Hamiltonian 
H, the stochastically evolved p is still a function of H, and so commutes with H at all times. 
Thus, for the choice of initial condition po = f{H), Eq. (10) simphfies to 

dp^]^a[{p,H}-2pTYpH]dWt . (11) 

Equation (11) defines p to be a martingale, for which the expectation Es conditional on 
information available up to time s obeys Es[pt\ = Ps , s < t, which reduces to E[p] = 
Eo[pt] — Po — f{H) when s = 0. [Note that if instead of Eq. (10) we had used the stochastic 
equation obtained from Eqs. (2) and (3b), the initial po = f{H) would not evolve in time at 
all, since Eq. (3b) vanishes identically when p commutes with H. This underscores again the 
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fact that Eqs. (3a) and (3b) are equivalent only for pure state density matrices, but define 
different stochastic evolutions for density matrices not obeying the pure state condition 
= P] 

Let us now show that at late times p evolves by Eq. (11) into a pure state projector 
when the Hamiltonian H is nondegenerate, or into a linear combination of projectors on a 
degenerate submanifold of H when H is degenerate. The proof of this parallels the proof 
[8,9] that Eqs. (la)-(ld) lead to state vector reduction. We consider the variance V of the 
Hamiltonian, defined by 

V = TrpH^ - {TrpHf , (12a) 

which by the Ito extension of the chain rule evolves in time as 

dV = TrdpH^ - 2TrpHTrdpH - {TrdpHf . (126) 

Using Eq. (11) for dp to evaluate TrdpH"', we find 

TrdpH"" = (7[Trp//"+^ - Tr pH^'Tr pH]dWt . (12c) 

Thus, substituting Eq. (12c) for n — 1, 2 into Eq. (12b) and taking the expectation, we get 

dE[V] = -a^E[V^]dt . (13) 

From here on the argument is identical to that of Refs. [8,9], and leads to the conclusion 
that as i — > oo the variance V approaches almost certainly. When the energy spectrum is 
nondegenerate, this implies that at late times only one density matrix element is nonzero, 
and so the initial density matrix po = f{H) has evolved to a pure state density matrix 
obeying = p. More generally, when the energy spectrum is degenerate, the vanishing of 
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the variance implies that the density matrix has evolved to a linear combination of projectors 
on a degenerate submanifold of the Hamiltonian. Thus, evolution of the initial density matrix 
po = f{H) by Eq. (10) leads to a late time density matrix that obeys E[p] — f{H), and 
which is a a pure state density matrix in the nondegenerate case, or a linear combination of 
projectors on a degenerate submanifold of H in the degenerate case. We take such density 
matrices as our model for the environment, and by the arguments of Sec. 2, are assured that 
the evolution of measurement systems uncoupled by Hamiltonian interaction terms to this 
environment are independent of the environmental dynamics, when the total system evolves 
under the density matrix dynamics of Eqs. (2) and (3a). 

4. Mean Field Approximation for a System 
Weakly Coupled to its Environment 

Let us now consider two subsystems with disjoint variables that are weakly coupled 
through an interaction term Aif in the Hamiltonian, so that the total Hamiltonian appearing 
in Eq. (10) is H = Hi + H2 + AH. We shall take subsystem 1 to be a measuring apparatus 
(including the microscopic system being measured), whose reduction dynamics we wish to 
follow, while we take subsystem 2 to be the external environment with which this measuring 
apparatus interacts. We shall derive a mean field approximation to the dynamics, in which 
each subsystem obeys an independent system stochastic equation with a modified Hamilto- 
nian, that reflects the mean interaction with the other subsystem. To this end, we substitute 
the independent subsystem Ansatz p — pip2 into Eq. (10), and take the partial trace Tr2 to 
average over the subsystem 2 dynamics, giving an effective equation for subsystem 1, and 
similarly, with the roles of 1 and 2 interchanged, to get an effective equation for subsystem 2. 
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We shall assume that in the limit of vanishing coupling Aif, the environment subsystem 2 
is in one of the ensembles constructed in Sec. 3 that is a function solely of H21 so that in the 
presence of Aif we have [p2,^^2] = O(Aif). We do not make a corresponding assumption 
for subsystem 1, since we will be interested in the case in which this is initially in a generic 
pure state. 

We proceed with this calculation term by term. From the left hand side of Eq. (10), 
substituting Eq. (6a) we get 

Tr2(ip = Tx2p2dpi + {pi + dpi)Tx2dp2 = dpi , (14) 

where we have used the condition Tr2P2 = 1 which implies that Tr2dp2 — 0. Prom the first 
term on the right hand side of Eq. (10) we get 

TT2{-i)[H, p]dt = -iTr2P2[Hi, pi] - iTr2[AH, pipa] - iTr2[H2, p2]pi ■ (15a) 

The first term on the right of Eq. (15a) gives simply 

-i[Hi,pi]dt . (156) 

Since Tr2Ai7p2 = Tr2P2Ai7, the second term on the right of Eq. (15a) becomes 

-i[Tr2P2AH,pi]dt , (15c) 

and the third term on the right of Eq. (15a) vanishes. So in sum, the first term on the right 
hand side of Eq. (10) gives 

-i[Hi + Tr2P2AH,pi]dt . (16) 

We turn next to the second term on the right hand side of Eq. (10), which gives 
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— ^a^dt times the partial trace of the double commutator, 

Ty2[H, [H, p]] = Ty2[H, + H2 + AH, [H, + H2 + AH, p,p2]] 
=Tr2[i/i + AH, [H, + H2 + AH, p,p2]] 

=Tr2{[i/i + AH, [H,, p,]p2] + [H,, [H2, p2]pi]] + [H,, [AH, pips]] (17) 

+[AH, [H2, P2]pi]] + [AH, [AH, P1P2]]} 

^[H^ + TY2P2AH, [H, + TY2P2AH, pi]] + 0{{AHf) , 

where we have used the facts that (i) Tr2[i?2, fl'(l) 2)] = for any function g of variables 1,2, 
and that (ii) by our equilibrium assumption for the environment, [if2,P2] is of order AH. 
[Step (ii) is the only one which does not go through in the corresponding effective equation 
calculation for the environment subsystem 2, leading to an additional term in its effective 
equation of motion given in Eq. (20a) below.] 

Finally, we turn to the third term on the right hand side of Eq. (10), which gives 
\adWt times 

I¥2[{p, H) - 2pTrpH] = Tr2[{piP2, Hi + H2 + AH} - 2pip2TriTr2PiP2(/^i + H2 + AH)] 
=Tr2[{piP2, + AH} + pi{p2, H2} - 2pip2Tripi{Hi + Tr2P2Aif) - 2pip2Tr2P2^^2] 
={pi, H, + TrsPsAi/} - 2piTripi(i/i + TrsPsAiJ) , 

(18) 

where no approximations have been made. 

Putting everything together, we see that the mean field approximation for the "mea- 
surement" subsystem 1 is 

dpi = -z[H[,pi]dt-la'[H[,[H[,pi]]dt+la[{pi,H[}-2piTTipiH[]dWt + 0{a\AH)^dt) , 
o Z 

(19a) 
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with the effective Hamiltonian 

H[^Hi+Tr2P2AH . (196) 

The corresponding equation for the "environment" subsystem 2 is obtained by interchanging 
the labels 1 and 2, and restoring the term dropped in step (ii) leading to Eq. (17), giving 

dp2 = - i[H'„ p2]dt - \a'[H'„ [H^, p2]]dt + la[{p2, H',} - 2p2Tr2P2mdWt 

(20a) 

--a'[Tr,{AH[H[, p,]), p^jdt + 0{a\AH))'dt) , 
with the effective Hamiltonian 

H^^H2 + TnpiAH . (206) 

The added term on the second hne of Eq. (20a) vanishes through order (Aif)^ when the 
reduction process for subsystem 1 has concluded, since then the density matrix for subsystem 
1 obeys = up to error terms of order (Aif)^. As a consistency check on the 

calculation, we see that the mean field evolution equations obey Tridpi = TY2dp2 — 0, and 
so preserve the trace conditions TVipi = Tr2P2 = 1- 

5. Dynamics of the Measurement Process: 
When is Coherence Maintained, When Does the State Vector Reduce? 

Let us now examine the implications of Eqs. (1-3) for measurements. We first have 

to specify the value of the parameter a governing the magnitude of the stochastic process. 

If quantum mechanics is modified at all, it seems likely that such modifications come from 

new physics at the Planck scale, and so we adopt for this discussion the estimate [7, 12] 
_ 1 

(7 ~ -^pi^ck' with -^Planck the Planck mass (in units with h — c — 1). With this estimate, 

15 



the reduction time tR in seconds for a state with initial energy variance AE is given [7,8,9] 



Thus, for AE equal to a proton mass, tji ~ 10 ^sec, while for AE equal to the mass of a 
nitrogen molecule, one has ~ 10~®sec. 



In order for stochastic energy-driven state vector reduction to give a viable phe- 
nomenology, it must satisfy the twin constraints of predicting the maintenance of coher- 
ence when this is observed, while predicting a rapid enough state vector reduction when 
a probabilistic choice between alternative outcomes is observed. We first discuss the con- 
straints imposed by the maintenance of coherence. We begin by noting that according to 
Eq. (21), the sole criterion governing how rapidly the state vector reduces is the energy 
variance; whether the system is microscopic or macroscopic plays no role. Coherent super- 
positions of macroscopic states, involving large numbers of particles, will persist in time if 
the energy spread between the superimposed states is small enough. As a first example, 
consider the recent superconducting quantum interference device (SQUID) experiments [13, 
14] observing the existence of coherent superpositions of macroscopic states consisting of 
oppositely circulating supercurrents. Taking for discussion the experiment [13] (which of 
the two has the larger energy variance between the superimposed states), the energy spread 
AE is roughly 8.6 x 10~^eV, and the circulating currents each correspond to the collective 
motion of ~ 10^ Cooper pairs. According to Eq. (21), despite the macroscopic structure 
of the state vector, the state vector reduction time tR for this experiment should be about 



by 




(21) 



5A. Maintenance of Coherence 
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10^^ sec ~ 3 X 10^^ years, and so maintenance of coherence is expected. 

As our next example of the maintenance of coherence in macroscopic systems, we 
consider a recent experiment [15] demonstrating diffraction of the fuUerenes Cqq and C70. 
We begin by noting that a diffraction pattern can be observed in a monoenergetic beam (in 
fact, this is the ideal condition for the experiment), so this class of experiments provides no 
evidence for coherent superpositions of states of differing energies. However, in a realistic 
experiment there will be an energy spread in the wave packet for each particle constituting 
the beam. To see a diffraction pattern, the spread in de Broglie wavelengths AA should be 
considerably smaller than A; adopting the very weak bound A A < A, we get the requirement 
that the spread in beam momenta Ap in each wave packet should obey Ap < p. This implies 
that each wave packet must have an energy spread AE obeying AE < 2£'kinetic- In the 
experiments of Ref. [15] the beam was obtained from an oven at approximately 900i^, and 
so the the bound on the energy spread becomes AE < 2 x (3/2) x 900K ~ .23eV. The 
corresponding state vector reduction time predicted by Eq. (21) is of order 1.5 x lO^'' sec ~ 
5 X 10^ years, and so energy driven state vector reduction plays no role in this experiment. 
Similar estimates, and the same conclusion, would hold if larger objects, such as viruses, 
were diffracted. 

These estimates suggest that in order to try to see the breakdown of coherence 
predicted by Eq. (21), one should consider experiments with systems having long lived 
metastable states separated by a large energy gap from the ground state. In atomic systems, 
the requirements on stability of the metastable state are very severe, since for a typical 
atomic energy splitting of a few eV, Eq. (21) predicts a state vector reduction time of order 
10^^ sec ~ 3 X 10^ years. For example, in the quantum intermittency experiments discussed 
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in [16,17], the metastable state lifetime is of order 1 sec, and so stochastic state vector 
reduction effects are neghgible. A potentially more promising case is provided by certain 
long-lived nuclear isomers [18], which are rendered metastable by their high spins, and which 
have large energy gaps from their ground states. For example, ^^^Hf has an isomer with a 
half life of 31 years suspended 2.4 MeV above its ground state. Quantum mechanics predicts 
that a coherent superposition of the isomeric state and the ground state should be stable for 
time intervals that are short relative to 31 years, whereas Eq. (21) predicts a spontaneous 
reduction of such a superposition to either the isomeric state or the ground state, with a 
reduction time of order 1 sec. The only nuclear isomer to exist naturally on earth is the 
metastable isomer of ^^'^Ta, which has a half-life of more than 10^^ years, an energy gap 
of 75keV from the ground state, and which accounts for roughly 1 part in 10"^ of naturally 
occurring tantalum. According to Eq. (21), a coherent superposition of the ground state and 
metastable isomer of ^*°Ta should spontaneously reduce to either the isomeric state or the 
ground state, with a reduction time of order 23 minutes. Maintenance of coherence of such 
a superposition for times significantly longer than this would decisively rule out Eqs. (1-3) 
as a phenomenology for state vector reduction. For example, if a laser using isomeric ^^°Ta 
could be constructed, and if the characteristic relaxation times for conventional sources of 
dissipation could be made much longer than 23 minutes, then the effects of Eq. (21) might 
appear as an additional, unconventional source of stochastic fluctuations or of dissipation. 
It would clearly be of interest to work out the detailed implications of Eqs. (1-3) for laser 
action in such a system. 

5B. Reduction in Measurement Situations 

We turn now to the second requirement that must be satisfied by a phenomenology 
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of state vector reduction, which is that it should lead to rapid reduction in experimental 
situations where a probabilistic outcome is observed. According to the von Neumann model 
for measurement [19], a measurement sets up a correlation between states 1/^) of a quantum 
system being measured, and macroscopically distinguishable states of the measuring 

apparatus Al, in such a way that an initial state 



An objective state vector reduction model must then account for the selection of one of the 
alternatives \ft)\M.t) from this superposition, with a probability given by |Qp. If the energy 
spread among the states |/^) is a typical atomic magnitude of a few eV, then as we have seen 
above using Eq. (21), the energy driven model of Eqs. (1-3) cannot quantitatively account 
for state vector reduction, unless the energy spreads among the alternative apparatus states 
in the superposition are much larger. Since in the ideal measurement model there is no 
energy transfer from the microscopic system to the apparatus, such an energy spread in the 
measurement apparatus states can only be present if induced by environmental interactions, 
which are ignored in the von Neumann analysis. If these environmentally induced energy 
fluctuations are large enough for the state vector to reduce in a time much smaller than the 
measurement time, then the observed results will agree with the Copenhagen interpretation 
of the measurement process. If the reduction time were of the order of or larger than the 
measurement time, then Eqs. (1-3) would predict stochastic fluctuations among alternative 
measurement outcomes, lasting until one is finally selected in a time roughly equal to Ir. 





evolves unitarily to 




(226) 
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However, as long as the apparatus states \Aie) are orthogonal for different £, no quantum 
interferences between different outcomes are possible. 

To reiterate, for environmental interactions to be effective in producing state vector 
reduction, they must lead to energy fluctuations of the apparatus in the course of a 
measurement, that are large enough for Eq. (21) to predict a reduction time that is less 
than the time it takes to make the measurement. Although different measuring devices have 
different response times, we shall for assume for purposes of our discussion that relevant 
measurement times range down to 10~^ seconds, which requires for reduction a /S.E ranging 
up to ~ 30GeV. We shall consider three possible sources of energy fluctuations: thermal 
energy fluctuations, fluctuations in apparatus mass from particle accretion processes, and 
fluctuations in apparatus mass from amplifled fluctuations in the currents that actuate the 
indicator devices. 

Thermal energy and temperature fluctuations in a canonical ensemble, that is, with 
flxed particle number, are governed by the equations 

{{i^Ef)Av^kBT''Cv , {{ATf)Av^kBTyCv , (23) 

with ks Boltzmann's constant and with Cy the heat capacity. Prom these formulas, and 
the values of the heat capacity and thermal conductivity for various substances, together 
with the formulas governing the surface radiation rate, one can estimate that when a body 
is large enough for the energy fluctuations at room temperature to be of order 1 to 30 GeV, 
the thermal relaxation time over which such energy fluctuations occur is much larger than 
measurement times of interest. (For example, for 1 gm of water at room temperature, the root 
mean square energy fluctuation is ~ 14GeV, and the thermal conduction relaxation time is > 
200 sec.) The reason for this is that the rate for heat transfer processes is proportional to the 
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temperature difference AT, and Eq. (23) shows that when a body is large, the temperature 
fluctuation AT is small. Thus we find that thermal energy fluctuations do not in general obey 
the criterion stated above for relevance to state vector reduction, that the energy fluctuation 
should occur within the measurement time. 

A more significant source of energy fluctuations comes [7] from particle accretion 
processes, for which we formulate a simple model within the mean field framework of Sec. 4. 
Consider a measuring apparatus which has N surface accretion sites for molecules of mass 
m. In Fock space representation, its Hamiltonian can be written as 

N 

Hi ^ Ho + 'ma'jaj , (24a) 

with Hq the bulk Hamiltonian for the apparatus, and with Oj and aj respectively the creation 
and annihilation operators for the accreted molecules. We assume the environment to contain 
a large number M of molecules that can be accreted onto the surface, with creation and 
annihilation operators fefc ,k — 1, M, and with a coupling to the accretion sites given 

by 

N M 

AH = Y[Ajka]bk + adjoint] . (246) 
j=i k=i 

This interaction Hamiltonian conserves the total number operator 

AT M 

7V = ^a]a,. + J]6l6, , (25) 

j=i k=i 

in other words, the total number of molecules accreted onto the surface or remaining in the 
environment is constant. 

In typical measurement situations, the environment density matrix will be diagonal 
in the number operator X^fli ^I^fe molecules being accreted. In this case, which we 
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term "incoherent" , the environmental expectation of AiJ vanishes, 

Tr2P2AH^0 , (26) 

and the reduction process is governed, according to Eqs. (19a,b), by the measurement system 
Hamiltonian Hi alone. (For a discussion of the coherent case, in which the environmental 
expectation of AH is nonzero, see the Appendix.) The Hamiltonian AH still plays a role, 
since in order AH in probability amplitudes [corresponding to order (Aif)^ in probabilities 
or transition rates] it leads to a sticking probability and an evaporation probability per unit 
time, respectively, for a molecule in the environment to accrete to the surface of the appa- 
ratus, and for a molecule already accreted to evaporate. As a result of these nonvanishing 
transition probabilities, the number of molecules accreted to the surface is constantly fluc- 
tuating. Assuming a simple colloid statistics model [20] in which each accretion site can 
hold only one molecule, the number of accreted molecules n obeys a Poisson distribution 
a{n,X) = e~^X"-/n\ with the mean X proportional to the sticking probability and in- 
versely proportional to the evaporation rate, and with the the root mean square fluctuation 
in the number of accreted molecules equal the square root of the mean number X of accreted 
molecules. 

Since distinguishable measurement outcomes must involve different configurations 
of the apparatus with respect to its environment, they will have different values of the 
accretion numbers Ojaj associated with the N accretion sites. Thus, the energy eigenvalue 
Hi of the measurement apparatus will differ for each distinguishable measurement outcome, 
with the spread of eigenvalues between any two outcomes being typically the mass of the 
accreted molecules m times the root mean square fluctuation in the number of accreted 
molecules. This statement assumes that the flux of accreting molecules in the environment 
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is high enough for such a fluctuation to actually occur during the state vector reduction 
time. In estimating when this condition holds, we will follow the review of Redhead [21] in 
assuming that the sticking probability is of order unity, in which case the minimum time for 
one molecule to be accreted onto an area of 1 cm^ can be read of from the molecular flux 
versus pressure tabulated in Table 2 of [21]. At room temperature and atmospheric pressure 
(760 Torr) the time for one molecule to be accreted onto an area of Icm^ is 3 x 10~^^ sec, 
while at an ultrahigh vacuum of 10~^^ Torr it is 3 x 10~^ sec. Thus, for an apparatus in the 
atmosphere at standard temperature and pressure, where the bulk of the accreting atoms are 
nitrogen molecules, the minimum apparatus area required for one molecule to accrete in a 
reduction time of 10~^sec (corresponding to a A£' equal to the mass of a nitrogen molecule) 
is 3 X 10~^^cm^, with the corresponding minimum area needed at a pressure of 10~^^ Torr 
equal to 3cm^. 

According to [21], the nighttime pressure at the surface of the moon is about 10~^^ 
Torr, while the pressure in interstellar space (within the galaxy) has been estimated as 10~^^ 
Torr. Under the assumption of a sticking probability of order unity, the mass accretion rate 
scales as the pressure divided by the mean molecular velocity. While molecular velocities 
away from the vicinity of the earth vary over a wide range, with effective temperatures in 
interstellar space ranging [22] from typically 50-100 to 10^ degrees Kelvin, we can get an 
estimate that is high by at most a factor of 2 or 3 by neglecting the velocity factor, and simply 
assuming that the mass accretion rate scales with the pressure from the values given in the 
table of [21]. With this assumption, the minimum apparatus area needed for a reduction 
time of 10~^sec is 3cm^ at the surface of the moon, and is 3 x lO^cm^ = 30m^ in interstellar 
space. In intergalactic space, the predominant matter [22] is highly ionized hydrogen, with 
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an effective temperature of order 10 degrees Kelvin and a density of ~ .23 proton per cubic 
meter. In this environment, the minimum apparatus area needed for a reduction time of 
10~^sec is around 8 x lO^m^ (corresponding to the accretion of 28 protons). If we were 
only to demand reduction in 3 x 10"'' second, then the needed apparatus area in interstellar 
space would decrease to less than lOcm^, while that in intergalactic space would decrease to 
~ Im^. Thus, a capsule large enough to sustain Schrodinger's cat, situated in intergalactic 
space, would have a reduction time stimulated by collisions with molecules in the intergalactic 
medium much smaller than the length of time needed to ascertain whether the cat were dead 
or alive! Perhaps more to the point, in a typical high precision molecular beam experiment 
[23], the beam velocity is of order lO^cm/sec, and the beam length is of order 2.7 meters. 
Hence the time for the beam to traverse the apparatus is 2.7 x 10~"^sec, and so the reduction 
time in intergalactic space for a capsule large enough to enclose the apparatus would be 
smaller, by at least an order of magnitude, than the measurement time. Clearly, in this 
situation the hmits predicted by Eq. (21) are being pushed, and there could be realizable 
experiments which, in intergalactic space, would be predicted to start to show evidence of 
the stochastic fluctuation between outcomes characteristic of the time evolution of the state 
vector in stochastic reduction models. But it seems unlikely that such an experiment could 
be devised within the confines of the solar system - the ambient matter fluxes are too high. 

In making some of the above estimates, it is convenient to have an alternative form 
of Eq. (21) that takes into account the accretion rate hmit on A£', and which is derived as 
follows. Let A4 be the mass accretion rate per unit area of the apparatus in units sec^\ so 
that the mass accretion on area A in time tR sec is AE — AM.tR. Substituting this into 
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Eq. (21) and solving for tR gives 



( 



2.8MeVy 
AM ) 



(27) 



This formula can be used whenever at least one molecule is accreted in the time tR. Given 
A4, we can calculate the area A corresponding to a given reduction time, and vice versa. 
For example, from Eq. (27) we find that for an apparatus of area Icm^ in the atmosphere at 
standard temperature and pressure, the reduction time is tR — 5 x 10~^^sec, corresponding 
to the accretion of ~ 1.5 x 10^ molecules in the time tR. 

Throughout this analysis, we have assumed that the Hamiltonian that is relevant for 
the stochastic Schrodinger equation is the total Hamiltonian 



defined by gravitational couplings to the stress-energy tensor T'^j,(x), which includes rest 
mass terms. Although in non-relativistic quantum mechanics one often drops rest mass 
terms when they lead to irrelevant constant energy shifts, there is no reason in principle to 
do so. In fact, in the standard model of elementary particles, all fermion rest masses arise 
from the Yukawa couphngs of the fermions to the Higgs particle, so that from this point of 
view rest masses are not an additive constant term in the Hamiltonian, but are a dynamical 
product of interactions. We have also assumed that the relevant surface area is that of 
the whole apparatus, rather than just that of components of potentially small area such as 
solid state detectors, emulsions, or particle collector cups. This assumption is motivated 
by our decoupling analysis of Sec. 2, where we saw that only noninteracting systems can 
be assumed (under certain equilibrium conditions) to decouple. The components of an 
apparatus (power supplies, magnets, vacuum pumps, detectors, indicator pointers, magnetic 




(28) 
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recording domains) are not in equilibrium and are in interaction with one another, and so 
using the area of the whole apparatus, rather than of just the smallest components, seems 
justified. 

We turn finally to a third potential source of energy fluctuations, arising from the 
amplified fluctuations in the currents which actuate experimental indicating or recording 
devices. Of course, if power sources are included, there are no overall current fluctua- 
tions, but power supplies are typically large in area and so when included in the system 
the accretion analysis just given indicates rapid reduction times. In a typical electrically 
amplified measurement, a final total charge transfer Ne (with e the charge of an electron) 
actuates an indicator or recording device. Assuming that the fluctuation in the current is 
the amplified fluctuation in the initially detected signal, for amplification gain G we have 
AA?" ~ G X [N/ G) 3 = {NG) 3 , an estimate which agrees within factors of order unity with 
the standard noise estimate for photomultipliers [24]. Let us take N to correspond to a 
charge transfer of 1 milliampere (a voltage change of 10 volts at 10 kf2 impedance) over a 
10~^ second pulse , so that A" ~ 6 x 10^, and assume a gain G ~ 10^, giving A A" ~ 8 x 10^. 
Multiplying by the electron mass of .5 x 10~^GeV, we find that the corresponding energy 
fluctuation is A£' ~ 4 x lO^GeV, which leads to state vector reduction in 5 x 10~^^sec. 
Thus, electric current fluctuations play a significant role in state vector reduction when the 
"apparatus" is defined to exclude power sources. 

Our overall conclusion is that conditions under which laboratory experiments are 
performed, as well as conditions under which space capsule experiments might be performed 
in the foreseeable future, are consistent with state vector reduction times as estimated by 
Eq. (21) that are well within experimental measurement times. 
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5C. Experiments with Semi-Silvered Mirrors 



In the preceding two subsections we have considered the case in which the energy 
variance is so small that coherence is maintained, and the case in which the energy variance 
is large enough that reduction proceeds rapidly to completion. Let us now briefly consider a 
case that contains elements of both, in which an apparatus is constructed using semi-silvered 
mirrors (for photons) or thin detectors (for particles) , so that there is a probability amplitude 
a for no interaction with the apparatus and the maintenance of coherence, and a probability 
amplitude j3 for a measurement to take place. The state vector then has the form, after the 
measurement interaction. 



Assuming that initial) and \Aie) differ sufficiently in energy for reduction to take place, 
there are now two possible classes of outcomes. With probability |/?Qp the £th measurement 
outcome is observed, while with probability jap the initial state is unchanged, corresponding 
to transmission through the semi-silvered mirror or thin detector. In making the latter 
assertion, we are using the fact that, as shown in Ref. [9], the model of Eqs. (1-3) obeys 
the Liiders projection postulate. That is, a component of the wave function lying within a 
submanifold of Hilbert space that is energy degenerate (or nearly degenerate, in the sense of 
Sec. 5A) survives unchanged in form, with the appropriate probability, as an outcome of the 
reduction process. This corresponds exactly to what happens when a beam is transmitted 
through a partially silvered mirror or a thin detector. This discussion generalizes immediately 
to the case in which the transmitted beam has different phase shifts in the various terms in 
the superposition over wave function components |/^). 




(29) 
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6. Coexistence of Reduction and Decoherence 

In the previous section, we considered the effects on an apparatus of inelastic coUi- 
sions, in which the mass of the apparatus fluctuates. However, an apparatus considerably 
more frequently suffers elastic collisions with atoms and photons in its environment, which 
are responsible for the decoherence effects [25,26,27] that have been much discussed in the 
literature. We shall argue in this section that decoherence effects do not substantially modify 
the results of the preceding section. 

We begin with a general argument that is independent of the details of modeling 
decoherence. When elastic interactions with the environment are taken into account, the 
effective apparatus wave function has to be extended to include the wave functions of all 
particles with which it interacts during the reduction time tR. Since this extension of the 
apparatus definition increases its area, the rate of mass fluctuations is increased, and the 
effective reduction time estimated in the preceding section is, if anything, decreased. (Be- 
cause tu is in general shorter than the time for molecules or photons in the environment 
to collide with one another, we do not have to continue this enlargement of the apparatus 
another step to include the particles with which the decohering particles interact - they can 
be regarded as effectively noninteracting. For example, we estimated above that in air at 
standard temperature and pressure, the reduction time tn for an apparatus of area Icm^ is 
~ 5 X 10~^^sec (during which time it accretes ~ 1.5 x 10^ air molecules), whereas the mean 
time between coUisions of air molecules with each other is ~ 10~"^°sec. Since the reduction 
time through accretion scales as the inverse | power of the density of the environmental 
medium, while the time between collisions of a molecule scales inversely as the density, in 
more dilute environments this inequality gets stronger.) With the definition of the apparatus 
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wave function extended in this way, Eq. (22b) is modified to read 



(30) 



with l^'e^) the environmental wave function associated with the ^th apparatus state. We can 
now apply the analysis developed above, using Eqs. (1-3) to describe the stochastic reduction 
of the wave function of Eq. (30), with the conclusions reached above unaltered. The wave 
function of the extended apparatus remains unit normalized, and its density matrix remains 
a pure state density matrix, which stochastically evolves to an energy eigenstate with the 
rate given by Eq. (21). 

Decoherence effects manifest themselves by the exponential decay with time of the 
inner product (^^e^l^e^')' ^' ^ so that the environmental states associated with different 

measurement outcomes become rapidly orthogonal. Correspondingly, in the reduced density 
matrix for the original unextended apparatus, obtained by tracing out the environmental 
degrees of freedom, there is an exponential decay of the off diagonal matrix elements. For 
the superposition of energy eigenstates that is relevant for our discussion, the relevant decay 
or decoherence rate of the off diagonal reduced matrix element is given by [25,26] 



Here A^scatt is the number of scatterings by environmental particles in unit time, (...) is an 
expectation in the state of the scattering particle, and and S^i are the scattering matrices 
acting on this particle when it scatters on the respective components of the apparatus wave 
function with state labels I and I' . In our context, these apparatus states differ only by the 
addition of some number of accreted molecules, and so in a weak scattering approximation 
the product S\S^' can be approximated as Sm, with Sm the scattering matrix for an envi- 



D = 7V,eattRe[l - {S\St,)] 



(31a) 
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ronmental particle to scatter from the accreted molecules. In this approximation, Eq. (31) 
simplifies to 

= A^scattRe[l - (5m)] . (316) 

Expressing the scattering matrix Sm in terms of the corresponding scattering amplitude, 
and using the optical theorem, Eq. (31b) reduces [25] to 

D = ^ X Scattering Rate . (32) 

Combining Eq. (32) with our estimates above for an area Icm^ apparatus in air, the deco- 
herence rate arising from molecules of the environment scattering on the molecules accreted 
during the reduction time is | x 10^°sec~^ x 1.5 x 10^~ .7 x lO^^sec"^, several orders of 
magnitude smaller than the reduction rate ~ 2 x lO^^sec"^, and so decoherence effects 
are in fact unimportant over the duration of the reduction process. Since the ratio of the 
reduction rate to the decoherence rate (calculated for the number of molecules accreted over 
the reduction time) scales with area as .45, this conclusion remains true down to an appa- 
ratus area of ~ 4 x 10~^^cm^, where the two rates become approximately equal. Thus, even 
for a very small apparatus the environmental states |^e^) are nearly identical to the initial 
environmental state |^e)) and so the wave function of Eq. (30) is negligibly entangled with 
the environment. Therefore instead of using the extended apparatus to discuss the reduction 
process, one is justified in ignoring decoherence and using the original unextended apparatus 
as in Eq. (22b). Our conclusion here for energy driven reduction models differs significantly 
from that reached [28,29] for spontaneous localization models, where decoherence effects over 
the reduction time are substantial; however, in these models the various apparatus states in 
the superposition of Eq. (22b) differ by a displacement of the center of mass of some part of 
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the apparatus, which gives a decoherence effect proportional to the (macroscopic) scattering 
cross section of that part of the apparatus which is displaced. 

7. Discussion and Conclusions 

The analysis we have given of a number of aspects of the effect of the environment on 
the measurement process, including the decoupling of isolated systems from environments 
in equilibrium, the effect of energy fluctuations induced by mass accretion, and the effect 
of decoherence processes, supports the view that the energy-driven stochastic Schrodinger 
equation gives a viable phenomenology of state vector reduction. According to this picture, 
a measurement takes place when the different outcomes are characterized by sufficiently 
large environmentally induced energy fluctuations in the apparatus for the state vector re- 
duction process, which is driven by the energy variance, to proceed rapidly to completion. 
The infinite Von Neumann regression (of an apparatus measuring an apparatus measuring 
an apparatus..., ad infinitum) terminates, when the apparatus size is large enough for its 
energy fluctuations to lead to state vector reduction within the specified observation time. 
This requirement on apparatus size meshes in a natural way with the intuitively obvious re- 
quirement that in a measurement, different experimental outcomes must be macroscopically 
distinguishable. 

Acknowledgments 

I wish to thank Jeeva Anandan, John Bahcall, Angelo Bassi, Todd Brun, Sudip 
Chakravarty, GianCarlo Ghirardi, Siyuan Han, Larry Horwitz, Lane Hughston, James Lukens, 
Indrajit Mitra, Ian Percival, Leo Stodolsky, and Prank Wilczek for stimulating discussions or 

31 



email correspondence, and to acknowledge the hospitality of the Aspen Center for Physics, 
where part of this work was done. This work was supported in part by the Department of 
Energy under Grant #DE-FG02-90ER40542. 



32 



Appendix. Coherent Case of the Accretion Model 

Anandan [30] has raised the interesting question of whether state vector reduction 
can proceed to a coherent state endpoint. One way in which this can happen is when the 
signal amplification process involves coherent states, as discussed in [27]. Another way, which 
we shall discuss here, corresponds to the "coherent" case of the accretion model formulated 
in Sec. 5B, in which the environment is in a coherent state, so that the environmental 
expectation of Aif is nonzero. Assuming for simplicity that there is only one accretion site, 
which can be multiply occupied, we have then 

Hi + T12P2AH ^ Ho + ma\ai + Xa\ + X*ai , (Ala) 

with A given by 

M 

A = ^AifcIV2P2&ik . {Alb) 

k=l 

Assuming Hq to commute with ai, Eqs. (Ala) and (Alb) describe the zero forcing 
frequency limit of the forced harmonic oscillator, which has been extensively studied [31,33], 
and can be succinctly solved by coherent state methods [32,33]. Defining z and Ci by 

z = —X/m , ci = ai — z = ai + X/m , (^420) 

we have 

Hi = Hq + mc\ci + constant , {A2b) 

which in its ci dependence is a standard harmonic oscillator. The ci oscillator ground state 
|0) obeys 

ci|0) = 0=^ai|0) =z|0) , {A2c) 
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in other words, |0) is a coherent state in terms of the original operators oi. 

Ignoring an overall constant arising from terms in Eq. (A2b) that commute with ai, 
the general eigenstate of Eq. (A2b) is \n), with n the number of ci quanta, and has energy 
eigenvalue mn. This state is a coherent superposition of states with different numbers of 
molecules on the accretion site. For energy eigenvalue n, the probability P{n\k) of finding 
n — k molecules on the site can be exactly expressed [31,33] as a Laguerre polynomial, and 
for \z\ « 1 and n large can be approximated [34] as 

P(n|A;)~[J|fc|(2ni|^|)]2 , {A3a) 

with Jfe the order k Bessel function; the Bessel function addition formula [35] 

oo 

1 = JoH' + 2^J„H2 {A3b) 

n=l 

implies that the probabilities of Eq. (A3a) sum to unity, 

oo 

P{n\k) = 1 . (A3c) 

k=—oo 

Equation (A3a) is rapidly oscillating as a function of k, but using the asymptotic estimate 
[36] 

/ 2 \ ^ 1 

Ji^iv sec 3) ~ cosiiytand — vd — -tt) , (^440) 

VTTi/tan/?/ V 4 / ' 

it is easily seen that the averaged envelope of P{n\k) is given by 

1 



(4n|2;p - k^)2 

showing that the values of k are mainly distributed (apart from an exponentially decaying 
tail) between —2n^\z\ and 2n2|2;|. 
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